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Parametric simulations with a Large Eddy Simulation (LES) model are used to explore 
the influence of crosswind shear on aircraft wake vortices. Previous studies based on field 
measurements, laboratory experiments, as well as LES, have shown that the vertical 
gradient of crosswind shear, Le. the second vertical derivative of the environmental 
crosswind ( d 2 ujdz \ can influence wake vortex transport. The presence of nonlinear 
vertical shear of the crosswind velocity can reduce the descent rate, causing a wake vortex 
pair to tilt and change in its lateral separation. The LES parametric studies confirm that the 
vertical gradient of crosswind shear does influence vortex trajectories. The parametric 
results also show that vortex decay from the effects of shear are complex since the crosswind 
shear ( dujdz), along with the vertical gradient of crosswind shear, can affect whether the 
lateral separation between wake vortices is increased or decreased. If the separation is 
decreased, the vortex linking time is decreased, and a more rapid decay of wake vortex 
circulation occurs. If the separation is increased, the time to link is increased, and at least 
one of the vortices of the vortex pair may have a longer life time than in the case without 
shear. In some cases, the wake vortices may never link. 
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Nomenclature 

wingspan of generating aircraft 

vortex separation 

initial vortex separation ( = nB/4 ) 

Brunt- Vaisala frequency 

crosswind shear gradient ( = d 2 ujdz 2 ) 

nondimensional S ( = S b 2 IV 0 ) 

time coordinate ( t = 0 at vortex generation ) 

nondimensional time ( = t V Q / b Q ) 

ambient crosswind 

initial vortex descent velocity ( = T 0 / (2 n b 0 ) ) 

vertical coordinate 

generating height 

altitude of port vortex 

altitude of starboard vortex 

vortex circulation 

initial vortex circulation 

nondimensional vortex circulation ( = T/T 0 ) 

eddy (turbulence) dissipation rate 

nondimensional eddy dissipation rate ( = (eZ? 0 ) 1/3 VJ l ) 

vertical shear of ambient crosswind 

vorticity of crosswind ( = -dujdz) 


I. Introduction 

Away from the influence of the ground, a wake-vortex pair typically descends due to its mutual induction while 
being transported horizontally with the environmental wind. The descent speed of the vortex pair decreases with 
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time as the circulation weakens or as the vortices become further apart. However, the presence of an environment 
with a vertical gradient in cross wind shear can alter the transport and decay of aircraft wake vortices. From field 
studies, 1 ’ 2 laboratory experiments, 3 ’ 4 and numerical simulations, 5 ’ 6 ’ 7 ’ 8 ’ 9 ’ 10 ’ 11 the vertical gradient of crosswind shear 
can preferentially reduce the descent speed of vortices, causing them to tilt, increase separation, and sometimes even 
to rise upward. Previous numerical investigations with the Terminal Area Simulation System (TASS) have shown 
this effect to be related to the second derivative of ambient crosswind, rather than the vertical shear of the 
crosswind. 8 This was also confirmed in the laboratory experiments conducted by Delisi and Robins. 4 

The schematic in Figure 1 illustrates how the tilt of a vortex pair is affected by the sign of the second derivative 
of the ambient crosswind. When the second derivative of the crosswind shear is negative, the counter-clockwise 
rotating vortex (depicted by a plus sign) rises relative to the position of clockwise rotating vortex (depicted by minus 
sign). The opposite occurs if the sign of the crosswind shear gradient is reversed. In this example, if the crosswind 
also increases with altitude, then the slower descending vortex would translate more quickly in the lateral direction 
than the faster descending vortex. This would cause an increasing vertical separation between the counter-rotating 
vortices and change the characteristic length scales that are based on the initial vortex circulation. The magnitude of 
descent velocity for each of the vortices would further decrease as the vortices move farther apart. In an opposite 
fashion, if the upwind vortex were to descend at a slower rate than the downwind as in the case of a positive 
crosswind shear gradient, the separation of the vortices would decrease. 



Figure 1. Schematic of crosswind shear effect on a wake vortex pair. Green curve represents vertical profile of 
crosswind when gradient of crosswind shear is either negative or positive. Assuming a right-handed coordinate 
system directed along the aircraft path, the port vortex is represented by the red circle, and the starboard by the 
blue circle. 


Although many of the effects of crosswind shear gradient on vortex transport are understood, little is known 
about how shear affects circulation decay. The observation of unequal vortex decay rates is not unusual during field 
measurement campaigns; and in some instances, the differential decay has led to the occurrence of solitary vortices 
that have prolonged lifetimes (e.g., B-757 Case 9 of the Idaho Falls test flights 12 ). 

Fast- time wake vortex prediction models are being advanced for operational applications as well as for use in 
safety-benefits studies. 13 Understanding of how crosswind shear affects wake vortex transport and decay is needed 
in order to improve the accuracy of the fast-time models. Most fast- time models ignore effects from the gradient of 
crosswind shear and assume both members of a vortex pair descend and decay at the same rate. However, in some 
instances the shear effects can be very important, since during unfavorable shear conditions wake vortices could 
linger within the flight corridor and persist over a longer period of time. 

In this paper we will use a numerical model to examine the behavior of wake vortices subject to idealized 
profiles of crosswind. In the following sections the setup and results of the parametric study are described. 
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II. Approach 

The Terminal Area Simulation System (TASS) is a three-dimensional large eddy simulation (LES) model that 
has been used in the past to study atmospheric and ground influences on wake vortex behavior. The TASS model is 
described in Refs [14, 15, 16] and validation studies are described in Refs [17, 18, 19]. The sensitivity of crosswind 
on vortex descent is examined from two sets of parametric experiments with TASS. Initial conditions are chosen to 
represent the wake vortices generated by a large category aircraft (B-737) and a heavy category aircraft (B-767). 
Environmental conditions are set to minimize the effects of turbulence (e* ~ 0.1) and stratification (N= 0) on vortex 
decay. Also, the wake vortices are initiated far enough from the ground to minimize any effect the ground surface 
may have. 

Sensitivity of vortex transport and decay is examined over a range of values for the gradient of crosswind shear. 
The vertical second derivative of ambient crosswind, or gradient of crosswind shear is defined by S = d 2 u c / dz 2 , and 
is assumed to be uniform within the initial domain. This parameter is nondimensionalized by the initial vortex 
separation, b Q , and the initial vortex descent rate, V Q , as S* = S b 2 / V Q . Experiments with two aircraft types are 
conducted since one of the aims is to develop a nondimensional relationship that is independent of aircraft type. 
Values assumed for S* in this study range from 0 to ± 0.5. Obviously if S* is a nontrivial constant, the value for 
dujdz will change with altitude. Therefore, experiments also are conducted for both dujdz > 0, and dujdz < 0. 

Turbulence for these experiments was grown by surface heating prior to adding the environmental crosswind. 
Tests were conducted to confirm that the initial crosswind profiles were not modified by the light turbulence during 
a three minute span equivalent to the lifetime of the wake vortices. The domain size assumed for the simulations 
was 450m x 650m x 455m in the longitudinal, lateral, and vertical directions, respectively. The grid size was set to 
1.5m in the longitudinal direction and 1.25m in the lateral and vertical directions. The lateral separation of the initial 
vortex pair, b Q , was assumed as 23m for the conditions representing the B-737 and 36m for the conditions 
representing the B-767. The initial circulation was assumed to be 250 m 2 s' 1 for the B-737 cases and 360 rns 1 for 
the B-767 cases. 


III. Results 

Results from the parametric simulations confirm that the crosswind shear gradient controls vortex tilt, with the 
port vortex descending slower for S*> 0, and the starboard vortex descending more slowly when S*<0. This is true 
irrespective of the sign of the crosswind shear. However, the crosswind shear, dujdz , is important in determining 
the rate at which the port and starboard vortex spread. For example, if S* < 0, the starboard vortex will be tilted 
upward from the port vortex; and if dujdz >0, the starboard vortex will be subject to a relatively stronger crosswind 
and will increase its lateral separation from the port vortex. Given the same S* 9 but now if dujdz <0, the starboard 
vortex would be exposed to a weaker crosswind than the port vortex, and a decreasing separation would follow. The 
TASS results show that the longest lived vortex is determined by the sign of the crosswind shear. If the vorticity 
associated with the crosswind shear is of opposite sign to the vorticity of the wake vortex, its lifetime will be shorter 
than the member who has vorticity with same sign as the crosswind shear. For example, if dujdz >0, the port vortex 
will live longer than the starboard vortex. These general conditions are summarized in Table 1. 


Table 1. Effect of crosswind shear on vortex pair. 


S* or 

dhijdz 1 

dujdz 

Tilt 

Change in Lateral 
Spread 

Longest Lived 
Vortex 

<0 

<0 

z s >z p 

decreasing 

starboard 

<0 

>0 

Z s >Zp 

increasing 

port 

>0 

>0 

Z s <Zp 

decreasing 

port 

>0 

<0 

Z s <Zp 

increasing 

starboard 

= 0 

^0 

Z s = Zp 

constant 

same 


The decrease or increase in the lateral separation due to the crosswind derivatives also plays an important role in 
the wake vortex decay. If the vortex lateral separation increases, then the characteristic length scale changes. For 
no shear cases this length scale is based on b Q . However, larger values may be more appropriate when the lateral 
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separation increases due to the crosswind shear gradient. An increasing vortex separation also would mean that 
vortex linking due to Crow instability 20 can be delayed or suppressed. In the case of decreasing vortex separation 
(Table 1), vortex linking would occur sooner. The nondimensional times when vortex linking occurred in our 
parametric study are shown in Tables 2 and 3. For conditions when the lateral spread increases, the linking times 
increase as well (Table 2). In the stronger shear cases, the separation becomes so great that linking does not occur. 
On the other hand, if the conditions are such that the lateral separation is decreased (Table 3), linking occurs sooner. 
It should be noted that dujdz and S*, affect the linking time since both control the rate at which the port and 
starboard vortices spread with time. 


Table 2. Effect of negative S* on nondimensional linking time when dujdz > 0. 


s* 

B-737 

B-767 

0.0 

3.0 

3.25 

-0.15 

3.0 

3.4 

-0.2 

3.25 

3.5 

-0.25 

3.35 

3.55 

-0.3 

GO 

GO 

-0.4 

GO 

GO 


Normally, the only environmental parameter thought to be important for predicting Crow linking is the 

01 90 O'! 

turbulence intensity, i.e., e*. ’ ’ Our results show that the presence of crosswind shear gradient can introduce great 
uncertainty into the prediction by either reducing or increasing the time at which of the vortex linking occurs. The 
ability to quantify the linking time is important in many wake vortex decay models. 


Table 3. Effect of positive S* on nondimensional linking time when dujdz > 0. 


S * 

B-737 

B-767 

0.0 

3.0 

3.25 

0.15 

2.1 

2.75 

0.2 

2.0 

2.45 

0.25 

1.9 

2.39 

0.3 

1.5 

2.1 

0.4 

1.4 

1.8 


The three-dimensional visual representation of the tilting effect due to the gradient of crosswind shear is shown 
in Figure 2. The vortex pair is spreading due to the crosswind shear gradient and tilting upward toward the starboard 
vortex. 

A visualization of one the cases with decreasing separation and accelerated linking time is shown in Figure 3. At 
the time of plot, the starboard vortex has nearly dissipated into rings, leaving an intact port vortex. 

Trajectory and circulation plots for two contrasting cases are shown in Figures 4 and 5. For the case with 
S*= -0.3, and dujdz > 0, the vortices never link and the port vortex has the longest lifetime. Even at three minutes, 
the port vortex has decayed only to 70% of its initial value, and remains many times stronger than the vortex in the 
case with no crosswind. For the case with S*= 0.2, and dujdz > 0, the lateral separation contracts, and Crow linking 
occurs at t = 21s. Note that both vortices decay faster than the baseline and the port vortex is slightly more 
persistent. In both Figure 4 and 5, the vortex trajectories are altered by the crosswind, and the vortices do not 
descend to the same depth as for the cases with no crosswind. 
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Figure 2. Three-dimensional perspective of the wake vortex pair tilted by shear. Viewed from the port side 
for the B-767 case with S*=- 0.25 and dujdz > 0, at t= 114 s (T*=5). Crosswind is to the left. 



Figure 3. Three-dimensional perspective of the wake with accelerated linking due to shear. Viewed from the 
port side for the B-737 case with S* = 0.2 and dujdz > 0, at t= 45s (T*= 3.4). Crosswind is to the right. 
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Figure 4. Vortex trajectory (top), and vortex circulation vs. time (bottom) for the B-767 case with S*= -0.3, 
and dujdz > 0. The baseline case for B-767 with no crosswind is shown as dashed line. The port vortex is 
denoted by red and the starboard vortex by the blue line. 
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Figure 5. Vortex trajectory (top), and vortex circulation vs. time (bottom) for the B-737 case with S*= 0.2, 
dujdz > 0 (same case as used in Figure 3). The baseline case for B-737 with no crosswind is shown by the 
dashed lines. The port vortex is denoted by red and the starboard vortex by blue. 


IV. Summary 

A parametric study is conducted with LES to examine the effect of idealized crosswind profiles on the behavior 
of a wake vortex pair. The crosswind shear gradient causes vortices in a vortex pair to descend at different rates. 
Once the vortices become tilted, they are transported at different lateral speeds due to the changes in crosswind with 
altitude (e.g., the higher vortex may be exposed to stronger crosswind). Their separation may increase or decrease 
depending upon the sign of the crosswind shear. Changes in the vortex separation due to the shear further affect the 
vortex descent rate, as well as the linking time and the rate of decay. Furthermore, the port and starboard members 
of the vortex pair may decay at different rates, and the vortex with the same sign of vorticity as the crosswind shear 
tends to persist longer than the vortex with opposite vorticity to the shear. The results of this study suggest that the 
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probabilistic approaches for wake prediction need to account for the greater uncertainty when shear is present, and 
deterministic approaches that predict time to link and vortex decay need to modify the characteristic scale due to the 
shear effects. 
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